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ABSTRACT 


This  report  covers  progress  made  in  the  second  year  of  the  resear  ...i 
program.  The  first  part  of  the  report  is  devoted  to  a  generalized  analyt¬ 
ical  prediction  of  the  ground  impulse  that  can  be  obtained  from  a  blast 
wave,  detonation  wave,  and  an  idealized  fuel-air  explosion.  The  latter 
consists  of  blast  wave  behavior  for  radius,  r,  less  than  a  critical  radius, 
r*,  and  Chapman -Jouguet  detonation  for  r  >  r*.  In  all  cases  so  far, 
the  finite  diameter  of  the  cloud  with  the  attendant  shock  wave  propagation 
beyond  the  cloud  has  not  been  taken  into  account.  The  latter  part  of  this 
report  is  devoted  to  the  experimental  aspects.  Improvements  in  the 
facility  for  generating  cylindrical  shock  waves  and  detonation  waves  are 
described.  Controlled  experiments  on  cylindrical  blast  waves  with  the 
associated  data  reduction  techniques  are  discussed.  The  results  are  in¬ 
terpreted  to  yield  a  calibration  on  the  effective  energy  release  of  the 
initiating  charge  of  Detasheet.  Two  phase  cylindrical  detonation  experi¬ 
ments  were  also  conducted  using  a  highly  refined  fraction  of  kerosene. 

The  results  indicate  that  at  small  radius  blast  wave  behavior  predominated 
whereas  at  larger  radius  a  constant  velocity  detonation  was  realized  when 
the  initiation  energy  was  sufficiently  high.  The  experimentally  determined 
transition  radius  between  the  two  types  of  behavior  agreed  very  well  with 
theoretical  values.  Cylindrical  detonations  in  gaseous  MAPP-air  mixtures 
were  also  studied.  The  variation  in  threshold  energy  required  for  initia¬ 
tion  as  well  as  rich  and  lean  limits  were  established.  The  results  agree 
very  well  with  large  scale  field  tests. 
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SECTION  I 


INTRODUCTION 

The  aim  of  this  research  program  is  to  gain  a  better  uaderetanding 
of  the  unconfined  explosion  of  detonatable  clouds.  Particular  attention 
is  given  to  the  liquid  fuel-air  explosion.  An  analytical  model  of  such 
explosions  which  would  be  capable  of  predicting  the  pressure  and  velocity 
field,  impulses  generated,  and  desirable  characteristics  of  the  fuel  has 
not  been  available.  Accordingly,  one  phase  of  this  research  has  been 
directed  to  gaining  such  an  understanding.  The  progress  on  this  phase 
is  described  in  the  next  section. 

On  the  experimental  side,  controlled  studies  of  cylindrical  and/  or 
spherical  heterogeneous  detonations  had  not  been  conducted.  Thus,the 
second  phase  of  this  research  is  concerned  with  experimental  studies 
on  a  laboratory  model  of  a  fuel-air  explosion.  This  work  is  described 


in  Section  in. 


SECTION  II 


THE  GROUND  IMPULSE  FROM  BLAST  WAVES,  DETONATIONS, 

AND  AN  IDEAL  FUEL-AIR  EXPLOSION 

A.  INTRODUCTION 

In  the  establishment  of  a  fuel-air  explosion  (FAE),  liquid  fuel  in  an 
appropriate  container  is  dispersed  into  the  atmosphere  as  a  cloud  of  fine 
droplets  by  the  detonation  of  a  primary  charge.  An  appropriately  placed 
and  timed  secondary  explosion  is  then  used  to  detonate  this  cloud  of  fuel 
droplets.  Methods  of  computing  the  ground  impulse  generated  by  an  ideal¬ 
ized  model  of  such  a  FAE  are  developed  in  this  report. 

The  shape  of  the  fuel  cloud  and  the  distribution  of  fuel  within  it 
which,  in  turn,  will  have  an  important  effect  on  the  ground  impulse,  will 
generally  be  quite  complex  and  depend  upon  the  method  of  fuel  dispersal. 
Here  an  idealized  FAE  with  cylindrical  symmetry  and  a  uniform  distri¬ 
bution  of  fuel  is  considered.  A  cylindrical  blast  wave  initiated  at  the 
center  of  the  cloud  sets  off  the  detonation.  This  ideal  ized  model  provides 
a  starting  point  for  the  consideration  of  other  FAE  phenomena. 

At  first  the  idealized  FAE  behaves  like  a  strong  cylindrical  blast 
wave  but  later  the  behavior  approaches  that  of  a  cylindrical  Chapman- 
Jouguet  (C-J)  detonation.  Universal  formulas  for  computing  the  total 
ground  impulse  generated  by  blast  waves  and  C-J  waves  are  first  developed 
below.  Then  it  is  show,  how  these  relations  can  be  combined  to  compute 
the  ground  impulse  generated  by  an  ideal  FAE.  Details  of  the  computations 
are  presented  in  Appendix  I. 
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B.  IDEALIZED  MODEL  OF  THE  FAE 


The  model  of  the  FAE  considered  here  consists  of  a  cloud  of  uni¬ 
formly  distributed  fuel  droplets  or  gaseous  fuel  adjacent  to  the  ground 
plane,  and  of  sufficient  height  so  that  the  effect  of  side  relief  can  be 
neglected.  The  secondary  blast  is  initiated  by  the  instantaneous  release 
of  explosive  along  the  axis  of  symmetry  of  the  FAE.  The  secondary 
blast  wave  and  the  subsequent  C-J  detonation  will  thus  be  cylindrical 
and  will  propagate  radially  outward  from  the  axis  of  symmetry.  The 
idealized  FAE  is  shown  in  Figure  1  and  is  equivalent  to  a  cylindrical.  FAE 
of  finite  height  but  confined  between  two  non-yielding  parallel  surfaces. 

It  has  been  possible  to  simulate  such  an  idealized  FAE  in  the  laboratory 
(Nicholls  et  alP). 

Blast  initiation  of  a  detonation,  which  corresponds  to  the  initial  phase 

/h\  /q\ 

of  a  FAE,  has  been  considered  by  KorobeinikoV'  7 and  Bach  et  al.  "'among 
others.  Initially  the  flow  is  dominated  by  the  strong  secondary  blast 
wave  while  the  energy  released  by  combustion  has  a  negligible  effect. 

The  blast  wave  decays  rapidly,  and  if  the  blast  energy  is  large  enough, 
a  C-J  detonation  is  established.  This  transition  from  blast  wave  to 
detonation  occurs  in  the  neighborhood  of  the  radius  r*  where  the  blast 
energy  E  is  equal  to  the  combustion  energy  release  within  r+,  i.  e.,  in 
the  region  r  <  r*.  The  radius  r+  plays  a  key  role  in  the  blast  initiation 
of  detonations  and  is  sometimes  referred  to  as  the  critical  blast  wave  radius. 
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Pressure  Distribution 


Figure  1.  Idealized  Fuel-Air  Explosion 


If  Q  is  the  combustion  energy  released  in  the  fuel  cloud  per  unit 


mass  of  mixture  and  if  Eq  is  the  secondary  blast  energy  released  per 
unit  length  on  the  axis  of  symmetry,  then  r*  will  be  defined  by  the  relation 

Eq  =  7r  r*2  Q  (1) 


where  is  the  initial  density  of  the  fuel  cloud.  Thus 

1/2 


r*  = 


Pi  Q 


(2) 


for  a  cylindrical  wave.  In  general,  with  v  =  1,  2,  and  3  for  plane,  cylin¬ 
drical,  and  spherical  symmetry 


r*  =  (v  EQ/av  Q  Pj) 


1/v 


(3) 


Eq  is  now  the  blast  energy  released  per  unit  area,  length,  or  the  total 
energy  released  at  a  point  in  the  spherical  v  =  3  case.  The  geometric 
factor  av  is  given  by 

=  2[v  -  1)  n  +  (v  -  2){v  -  3) 


and  a  =  2,  2v,  An  for  v  -  1,  2,  3. 

InitiaTy  when  r  r+,  where  r  is  the  radius  of  the  outward  propa- 
s  s 

gating  blast  or  detonation,  the  FAE  is  essentially  a  strong  blast  wave  and 

fa)  (5) 

can  be  described  by  the  self- similar  solution  of  Sedo^  'and  Taylor  . 

When  r  »  r*  the  FAE  becomes  a  C-J  wave  which  can  also  be  described 
s 


by  a  self-similar  solution  of  the  conservation  equations.  In  the  transi¬ 


tional  region  r  ~  r*  neither  the  blast  wave  or  the  C-J  similarity  solu- 
s 

tions  will  be  valid.  A  semi-empirical  theory  for  flow  in  this  transitional 

(3) 

region  has  been  developed  by  Bach  et  al : 

In  the  idealized,  model  of  the  FAE  the  complications  which  arise  in 

the  transitional  region  are  ignored  by  assuming  that  the  FAE  can  be 

described  by  the  self-similar  strong  blast  wave  solution  when  r  <  r+ 

s 

and  by  the  self-similar  solution  for  2  C-J  detonation  when  rg  >  r+.  This 
approximation  is  a  key  feature  of  the  idealized  FAE  and  leads  to  great 
simplification.  Experiments®  indicate  that  the  transition  from  blast  to 
detonative  behavior  occurs  very  rapidly  near  rg  ~  r+.  Thus,  neglect  of 
the  details  of  the  transition  should  result  in  only  minor  errors  in  the 
determination  of  FAE  propagation  and  impulse- 


It  can  be  seen  from  Equations  (2)  and  (3)  that  the  critical  blast  wave 


radius  r*  depends  on  both  the  secondary  blast  energy  EQ  and  the  heat 

release  Q  of  the  fuel.  Unless  E  exceeds  a  certain  minimum  initiation 

0 

energy  (Eo)cr-|.j  the  secondary  blast  will  decay  *.o  an  acoustic  wave  and 
the  fuel  cloud  will  fail  to  detonate.  As  indicated  in  Reference  1,  (E  )  ..  will 

depend  upon  the  structure  of  the  detonation  wave,  and  with  the  present 


state  of  knowledge  actual  values  of  (E0)cr„  can  only  be  determined  experi¬ 
mentally.  In  an  actual  FAE  device  Eq  >  (E0)cr^,  *n  general.  Precise 
determination  of  Q  requires  detailed  calculation  of  the  equilibrium 
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composition  of  the  combustion  products  behind  the  C-J  detonation. 

How  ever,  as  shown  in  Reference  1,  remarkably  accurate  results  can  be 
obtained  for  the  jump  condition  across  detonations  if  the  perfect  gas 
equation  is  used  both  upstream  and  downstream  of  the  C-J  discontinuity 
with,  however,  different  values  of  molecular  weight  and  the  ratio  of 
specific  heats.  Then 

C2 

Q  = - K -  (4) 

2(y22  -  1) 

where  C  is  the  velocity  of  propagation  and  y2  is  the  ratio  of  specific  heats 
of  the  combustion  products.  With  Equation  (4)  the  expression  for  r*  can  also 
be  written  in  the  form 


(5) 


The  self-similar  blast  wave  and  C-J  solutions  are  patched  together 
when  r  =  r*  in  the  idealized  FAE  model  described  above.  This  means 
that  the  entire  flow  field  r  <  rg  corresponds  to  that  of  a  blast  wave  when 
<  r+  and  to  that  of  a  C-J  detonation  when  x n  >  r*.  The  pressures 
P2  behind  the  blast  wave  and  p^j  behind  a  C-J  detonation  are  not  equal 
at  the  point  of  transition  r  =  r.  but  are  of  the  same  order  of  magnitude. 
From  the  simple  theory  for  C-J  detonations 
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The  parameter  a{y,  v)  is  a  function  of  y  and  v  and  is  plotted  in  Figure  2. 
Letting  p*  be  the  value  of  p^  corresponding  to  rg  =  r^  and  introducing 
the  expression  (3)  for  r^  in  Equation  (7)  then  yields 


4a 

_ v _ 

( v  +  2)va(y1  +  1  )(y2  -  1) 


(8) 


Taking  y ^  =  1.4.  y^  =  1*  2,  a  =  1.0,  Equation  (8)  yields  the  following  result 


V 

1 

2 

3 

p2*,/pCJ 

1.85 

1.64 

i,  39 

so  that  the  blast  and  C-J  pressures  are  indeed  of  the  same  order  of 
magnitude  when  rg  -  r+. 

In  order  to  determine  the  ground  impulse  generated  by  an  ideal  FAE, 
it  is  necessary  to  determine  the  ground  impulse  of  a  strong  blast  wave  and 
of  a  C-J  wave  as  described  below. 


8 


Figure  2.  Variation  of  the  Blast  Wave 
ot(y,  v)  with  y 


C.  TOTAL  GROUND  IMPULSE 


1.  Integral  for  the  Impulse 

The  pressure,  p,  within  the  leading  shock  of  a  FAE,  i.e-,  in  the  region 

r  <  r  will  be  a  function  of  the  radial  distance  r  from  the  blast  center 
s 

and  of  time  t.  If  p^  is  the  ambient  pressure  ahead  of  the  FAE,  then  the 
ground  impulse  generated  by  the  FAE  per  unit  area  at  a  given  point  dur¬ 
ing  time  interval  dt  is  (p  -  pj)  dt.  The  total  impulse  in  the  region  r  <rg 
generated  during  interval  dt  is  given  by  the  integral 


The  upper  limit  r  is,  as  indicated,  a  function  of  time.  Finally,  the  total 
s 

ground  impulse,  I{t),  at  time  t  after  the  initiation  of  the  secondary  blast 
on  the  axis  of  symmetry  will  be 


It  is  often  also  of  interest  to  determine  the  ground  impulse  generated 
by  a  planar  and  spherical  wave.  In  the  planar  case  the  secondary. blast 
is  initiated  on  the  plane  r  =0  where  now  r  is  the  distance  from  blast 
center.  In  the  planar  case,  v  =  1,  r  is  a  Cartesian  Coordinate  and  can 
have  positive  and  negative  values.  As  shown  in  Figure  3,  two  planar  blast 
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or  detonation  waves  then  propagate  symmetrically  in  both  the  r  and 


-  r  directions.  The  planar  configuration  Is  typical  of  that  encountered  in 
various  bag  tests.  The  total  ground  impulse  generated  by  a  planar  wave 


is  given  by 


CIO'* 


The  factor  of  2  in  Equation  (10)  accounts  for  the  propagation  in  the  -aid  -  r 
directions. 


The  pressure  generated  by  a  spherical  wave  with  the  secondary 


blast  initiated  a*,  a  point  on  the  ground  plane  (as  shown  in  Figure  4)  still 


acts  on  a  circular  area  so  that  the  impulse  is  still  given  by  Equation  (S). 
Now.,  however  p(r.  0  and  r  (J)  correspond  to  a  spherical  wave. 

In  general,  evaluation  of  u-e  impulse  from  Equations  (9) ana  (10)  will  re¬ 
quire  extensive  numerical  computation.  However,  in  the  case  of  strong 


blast  waves  and  C-J  detonations,  the  integrals  (3)  and  (10)  can  be  reduced 
to  simple  expressions  for  the  impulse  fit)  as  shown  below. 

2.  Ground  Impulse  for  a  Strong  Blast  Wave 

A  detailed  discuss. on  of  the  self-similar  solution  for  a  strong  blast 
wave  is  given  in  Reference  1.  For  a  strong  blast  wave  the  ambient  pressure 
«.  p,  and  so  can  be  neglected.  The  pressure  p  depends  on  the  dimen¬ 
sionless  radius  A  =  r  r  according  to 


Spherical  FAE — Quarter  View 


p  =  pi\7A(x,l/) 


P^(A,  v)  is  a  universal  function  of  A  which  depends  only  on  the  ratio  of 

specific  heats  >,  and  on  the  geometry  of  the  wave,  i.e.,  on  v.  The  radius 

r  (t)  of  the  shock  front  is  given  by 
s 


P1  a  (v,v) 


v+2  2 

,v+2 


v  =  1,  2,  3 


The  function  a{y,  u),  which  also  occurs  in  Equation  (7)  for  Pg,  is  shown 
ir  Figure  2. 

Neglecting  p^  and  introducing  the  variables  A  and  P  the  inner  inte¬ 
gral  in  Equations  (9)  and  (10)  becomes 


s  1  r  2  x2 

J  P  S'1  dr  =  <jf  PlJ-5~'Pb(KMrs\f~1rs 

A  A 


where  Equation  (13)  reduces  to  the  planar  and  cylindrical  forms  for  v  -  1,2. 
Upon  introducing  Equation  (12)  for  r ^(t),  Equation  (13)  becomes 


%/ 


/  p  r  dr  =  %J  Pb(X’  V)  X"+  dX 
0  0 
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The  integral  in  Equation  (14)  is  a  universal  function  of  y  and  v  and 
can  be  evaluated  from  the  blast  wave  solution,  as  shown  in  Appendix  I. 
Integration  with  respect  to  time  t  now  yields  the  follow  mg  result  for 
blast  wave  ground  impulse 


I(t)  =  a  E  t  ;  v  =  1,2 

'  v  oa(y,v)  *  ’ 


(15) 


p(y,  v)  =  f  J?^)u)xu+1  dk 


/ 


For  planar  and  cylindrical  blast  waves  the  ground  impulse  thus  varies 
linearly  with  time  and  is  proportional  to  the  blast  energy.  The  universal 
function  fi(y,  v)  is  plotted  in  Figure  5,  while  Figure  6  shows  the  varia¬ 
tion  of  a  J3/  a  which  is  equivalent  to  a  dimensionless  impulse  I(t)/ EQt. 

In  the  planar,  v  =  1,  case  T(t)  is  the  impulse  on  a  strip  of  ground  of 
unit  width.  The  values  of  a(y,v),  P(y,v),  and  cr ^3/ a  are  tabulated  in 
Table  I. 

For  a  spherical  blast  wave  which,  however,  acts  on  a  circular 
region  of  the  ground  as  shown  in  Figure  4,  the  results  are  somewhat 
different.  Equation  (11)  relation  p  to  P^,  v)  remains  valid.  Now, 
however,  the  inner  integral  of  Equation  (9)  becomes 


r  (t) 
s' ' 


2n 


f  pr  dr  =  277  /  p. 

r\  a  * 


2  2 
r  X 
s 


0 


P,  (X,  3)  r  X  r  dX 
b'  ’  '  s  s 


(16) 
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TABLE  I.  BLAST  WAVE 


(a)  Values  of a(y,  v) 

B 

v-  1.0 

i/  -  2.0 

BQ 

1. 1 

2.  9383 

2.6449 

2.2997 

1.2 

1.9114 

1.7622 

1.5383 

1.3 

1.362t 

1.2661 

1. 1005 

1.4 

1.0264 

0.  9666 

0.8428 

1.5 

0.8163 

0.7711 

0.  6743 

1.6 

0.6695 

0.6346 

0.  5572 

1.7 

0.5619 

0.5356 

0.4716 

(b)  Values  of  0(y,  v) 

y 

n 

o 

v  =  2.0 

v  =  3.0 

1.1 

0.20822 

0.05857 

0.  03669 

1.2 

0.19613 

0.05521 

0. 03397 

1.3 

0. 18557 

0.05230 

0.03169 

1.4 

0.17621 

0.  04972 

0. 02973 

1.5 

0. 16784 

0. 04743 

0. 02803 

1.6 

0. 16028 

0.  04536 

0. 02654 

1.7 

0. 154342 

0. 04349 

0. 02522 

(c)  Values  of  Dimensionless  Impulse  Function  vs  y 


a  m 

10  _  0(y) 

v  a  (y) 

3  "  .4/5 

a  (y) 

y 

v=l.O 

v  =  2.0 

v  =  3.0 

l.i 

0. 14173 

0. 13914 

0. 197348 

1.2 

0.20523 

0. 19687 

0.252606 

1.3 

0.27250 

0.25554 

0.  307340 

1.4 

0.34336 

0.32322 

0.356979 

1.5 

0.41121 

0.38544 

0. 402344 

1.6 

0.47884 

0.  44916 

0.443792 

1.7 

0.54611 

0.51020 

0.481910 
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Now,  however 


so  that  Equation  (16)  becomes 


(16a) 


After  integration  with  respect  to  time  the  spherical  blast  wave  ground 
impulse  becomes 


I(t)  =  ^  »px  (E0/oPl)4/b  (3(y,3)t3/5 


(IT) 


where  P(y,  3)  is  again  defined  by  Equation  (15),  but  nov  Pb(\,3)  corresponds  to 
the  spherical  blast  wave  solution.  The  function  £(y,  v)  for  v  =  3  and  the 
normalized  impulse 


I(t) 

Pl  <V»>1> 


10 

I7s  =  T1T 


P(r,  3) 


(18) 


are  plotted  in  Figures  5  and  6. 

Using  Equation  (12)  for  r  (t)  the  formulas  for  I(t)  may  also  be  expressed 

5 

in  terms  of  the  blast  radius  r  ,  a  form  more  useful  in  certain  applica- 

s 

tions.  Thus, for  v  =  1,  2 
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Knowing  a  and  8  from  Figures  2  and  5,  I(t)  and  I(r  )  can  now  be  readily 

*  b 


calculated  for  anv  E  and  p, . 

o  1 


3.  Ground  Impulse  from  a  Chapman -Jouguet  Detonation 

The  ground  impulse  generated  by  a  C-J  detonation  Is  also  given 
by  Equations  (9)  and  (10).  Now,  however,  the  ambient  pressure  p^  is  no 
longer  negligible  compared  to  the  pressure  behind  the  detonation  front. 

The  self  similar  solution  for  a  C-J  detonation,  which  is  discussed  in  detail 
in  Reference  1  can  be  used  to  reduce  the  impulse  integral  [  (Equations  (9) 
and(10))]  to  a  relatively  simple  form. 

As  before, a  dimensionless  pressure  P^X,  v)  is  defined  by 

2 

P  =  bi  ^2  pd  (x>  v)  (ID 

where  X  =  r/r 

s 

In  contrast  to  the  blast  wave,  a  detonation  front  propagates  with  a  constant 
velocity  C  so  that 

r  =  Ct  (21) 

s 

Introducing  the  similarity  variables  the  inner  integral  in  Equations  (9) 
and  (10)  becomes 
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]  (p  "  pl^  rI,_1  ^ 
0 


~r+2  .  y 
a  C  t 
y  1 


Pd  (A,  y)  A^+1  dA  - - y- 

Pi  C  y 


;  ^  -  l  >  2. 


(22) 


The  Mach  number,  of  the  detonation  is  given  by 


md  = 


^  (p l/pl) 


where  y^  is  the  ratio  of  specific  heats  of  the  unburned  fuel-  Hence,  the 
ground  impulse  I(t)  will  be  given  by 


I(t)  = 


„v+2  ,  y+1 
°v  P1  C 
y  +  1 


i. 

/ 

0 


Pd  (A,  y)  Ay+1  dA - - 


^1MD 


,  y  =  1,  2  (23) 


For  a  spherical  detonation  the  ground  impulse  is  given  by 


I(t)  = 


2tt  t3 


j  Pd  (A,  3)  A3  dA - - 


0 


2^i  md 


(24) 


The  integral  which  appears  in  Equations  (22)  and  (23),  and  which  is 
denoted  by 
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1 

5(y2, v)  -  J  pd(*> V)  xV+1  ^  v  =  1>  2 

o 

.  (25) 

1 

5(y2,3)  =  y  Pd(A,3)  A3  dA  v  =  3 

0 

depends  upon  the  ratio  of  specific  heats,  y^,  of  the  combustion  products, 
upon  the  geometry  of  the  detonation,  and,  weakly,  upon  the  detonation 
Mach  number,  «D-  The  computation  of  6(y2,  u)  is  discussed  in  detail 
in  Appendix  n. 

The  ground  impulse  also  can  be  expressed  as  a  function  of  the  det¬ 
onation  radius  r  so  that 
s 


a  p,  C  r 
v  1  s 

V  +  i 


i>+l 


5(y2,  y)  - 


’’’'l  “d" 


i/  =  1, 2 


and  for  spherical  waves 


2v  p1  C  r  1 

L  2yi  “d  J 


v=  3 


(26) 


(27) 


The  expressions  for  the  ground  impulse  generated  by  a  C-J  detonation 
are  summarized  in  Table  D3. 

The  detonation  velocity  C,  Mach  number,  M^,  ratio  of  specific 
heats  y2,  and  the  jump  in  pressure,  density,  and  temperature  across 
the  detonation  are  required  ju  order  to  compute  6(r2,  v)  and  the  ground 


23 


KrJ 


impulse  I(t)  or  I(r  )  from  the  above  relations.  In  other  words,  the 
s 

Chapman -Jouguet  conditions  must  be  computed  and,  of  course,  vary 
with  the  fuel  used  and  with  the  mixture  ratio.  Precise  determination  of 
the  C-J  conditions  is  tedious  since  the  equilibrium  composition  of  the 
_  combustion  products  must  be-  determined.  Fortunately  Gordon  and 

/C\ 

McBride'  *  have  developed  a  very  efficient  computer  program  fern  making 
such  calculations,  and  this  program  has  been  used  to  determine  the  C-J 
conditions  for  gaseous  air-MAPP  and  air -methane  mixtures  as  outlined 
in  Appendix  HI.  These  two  fuel -air  mi-lures  have  been  chosen  as  repre¬ 
sentative  of  the  hydrocarbon  fuel-air  mixtures  which  might  be  used  in  FAE's. 
The  detonation  parameters  for  methane  and  MAPP-air  mixtures  are  shown 
in  Figures  7  and  8  and  tabulated  in  Table  IV. 

The  function  6(y9,  v)  has  been  calculated  for  both  air-methane  and 
air-MAPP  mixtures.  For  the  two  fuels  considered,  ranges  between 
1. 16  and  1. 31  while  is  between  4. 0  and  5. 0,  and  over  this  range 
6(y2,  v)  is  almost  independent  of  y 2  and  the  mixture  ratio  #.  In  fact, 

6(72,  v)  depends  only  on  v  for  CH^  and  MAPP,  and  the  average  value  of 
6(y2,  v )  is  given  in  Table  V  below.  The  values  of  yn  and  are  typical 
of  most  hydrocarbon  air  C-J  detonations;  hence,  the  values  of  6  (y2,  v) 
in  Table  V  can  be  used  to  compute  the  impulse  for  hydrocarbon -air 
C-J  detonations,  in  general,  with  reasonable  accuracy. 
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Figure  8(a).  Detonation  Properties  —  MAPP-Air 


Figure  8(b) 


TABLE  IV(a).  PROPERTIES  OF  METHANE  AIR  DETONATIONS 


25.873  19.620  1.3088  1  10.70  4.80  4945.150  4.12  0.0020520 


TABLE  IV(b).  PROPERTIES  OF  MAPP-AER  DETONATIONS 

Air-Mapp  _ 


TABLE  V.  AVERAGE  VALUE  OF 


6(>2’  v) 


V 

6(y2-  ‘"'avg 

1 

0.  190 

2 

0. 138 

3 

0. 135 

The  quantities  for  v  =  1,  2  and  I(t)/t^  for  v  =  3  depend  upon 

Pp  C,  and  MD  [  (Equations  (23),  (24))]  and  so  vary  appreciably  with  the  equiv- 
lence  ratio  <f>.  With  these  impulse  functions,  which  are  plotted  in  Figures  9 
and  10,  I(t)  is  readily  determined  for  C-J  waves  in  air  methane  and  air- 
MAPP  mixtures.  It  is  of  interest  to  note  that  the  peak  values  of  I(t)  seem 
to  occur  for  rich  mixtures  with  an  equivalence  ratio  of  <f>  ^  1.  2  for  air 
methane  and  <f>  «  1.  5  for  air-MAPP. 

Since  r  =  Ct  it  follows  that 
s 


and 


(28) 


r 

s 


Mi 

t3  c3 


v  -  3 


(29) 
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Total  Ground  Impulse  CH^-Alr 
=  1  atm,  T1  =  298. 16°K 


Total  Ground  Impulse  Mapp-Air 
=  1  atm,  T,  =  298. 16°K 


4.  Ground  Impulse  from  Ideal  Fuel -Air  Explosion 

As  indicated  in  subsection  B,  above,  the  ideal  model  of  a  FAE  behaves 

as  a  blast  wave  for  r  <  r*  and  as  a  Chapman -Joug  ugt  Detonation  for 

s 

r  >  r+.  The  ground  impulse  generated  by  this  ideal  FAE  is  determined 
s 

below. 

A  wave  with  cylindrical  symmetry,  v  -  2,  is  considered.  Then  for 

r  <  r+  the  impulse  is  simply  that  due  to  a  blast  wave  and  is  gi\en  by 
s 

Equation  (19) 


When  rg  >  r+,  the  entire  FAE  flow  field  is  replaced  by  that  due  to  a 
Chapman-Jouguet  detonation.  Hence  the  ground  impulse  is  approximated 

fcy 

T<rs>  Vrs>  -TCJ(r.'  +TBW(r*'  (30) 

Thus  when  r  =  r+,  I(r  )  corresponds  to  the  blast  wave  ground  impulse, 
s  s 

while  with  r  >  r  ,  I(r  )  is  the  ground  impulse  generated  by  the  blast 
s  s 

wave  up  to  r  =  r+  plus  the  ground  impulse  generated  by  the  C-J  detona- 
s 

tion  for  rg  >  r+.  Subscripts  CJ  and  BW  refer  to  the  detonation  and  blast 
wave,  respectively. 

Using  the  relations  for  blast  wave  and  C-J  wave  impulse  Equation  (30) 
becomes 
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It  is  assumed  that  the  blast  wave  propagates  through  the  unburned  gas 
so  that  B(y,  2)  is  evaluated  for  y  -  Y^-  Using  Equation  (5)  the  expression  for 
r+  can  also  be  written  in  the  form 


2  2Eo  (y22  "  U 

r*  = - 2 

T  Pi  MD 


(32) 


The  results  above  can  also  be  expressed  in  terms  of  the  time  t  from 
the  initiation  of  the  secondary  blast.  Thus,  for 

t  <  t+  ;  I(t)  =  2ir  Eq  O/a)  t  (33) 


where  t+  corresponds  to  the  time  at  which  the  blast  wave  radius  reaches 
the  critical  value  r*.  Thus 


t 


* 


(34) 
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For  t  >  t+ 


I(t)  = 


2n  pj  C 


Hy2,  2)  - 


2^ 


(t3  -  o 


+  2n  Eq  (/3/a)  t, 


(35) 


5.  Sample  Computation 

At  this  point  it  is  useful  to  carryout  a  sample  computation.  Experi¬ 
ments  in  the  segmented  test  section  described  in  Reference  1,  indicated  that 
the  critical  value  of  Eo  to  initiate  a  C-J  detonation  in  an  air-MAPP  mix¬ 
ture  with  equivalence  ratio  <f>  =  0.  563  is 

Eq  =  4.38  X  105  (ft -lb) /ft 

which  is  equivalent  to  150  gm  of  Detasheet  per  foot  (using  an  energy  re¬ 
lease  of  2911  ft-lb  of  energy  release  per  gm  of  Detasheet^.  For  the 
above  air  -MAPP  mixture 

Md  =  4.  6S 
Y1  =  1.30 

yi  ~ 1'23 

Then  with  p^  =  14.  7  psia,  i.  e.,  the  sea  level  atmospheric  pressure 

5  I1/2 

r  =  2(4.  38)  x  10  lb  (1.51-1) 

3. 14(1.  30) (14.  7)(144) (lb/ft2) (21, 8) 

— 

=  1.54  ft  • 


P1  =  2.  29  x  1Q~3  slug/ft3 
C  =  5320  ft/sec 
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From  Figure  4  o(y1t  2)  =  1.  22  where  is  the  ratio  of  specific 
heats  of  the  unburned  fuel-air  mixture. 

From  Figure  5  |3(y,  v)  =  /3(1.  o,  2)  =  0.  052.  Then  with  r  <  r*  the 

blast  wave  impulse  from  Equation  (19)  becomes 


/E  p  \1/2 

I(rs> =  2l,(~ir)  3  r* 


.  6.28^- 


5  -3\1/2 

4.38  x  10  x  2. 29  x  10  \  2 

- L22 - )  (0-  0:>2>  rs 


=  9.  36  r  lb  sec 
s 


r  in  ft 
s 


At  the  critical  radius  r„ 


!BW(r  *)  =  9.  36  r+2  =  9.  36(1.  54) 2  =  22.  2  lb  sec 


The  next  step  is  to  compute  the  impulse  *Qj(rs)  generated  by  the 
detonation  through  the  air-MAPP  mixture.  I _  (r  )  is  given  by 

L-o  S 


TCJ(r  '  ■ 


2t7  p..  C  r 
1  s 


6('2’  21 


*1  V 


As  indicated  in  the  discussion  above, v)  is  essentially  independent 
of  >2  so  that  the  average  value  6 (y  2)  =  0. 138  can  be  used.  Then 


(26) 
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3 


s 

=  3. 06  r  3  lb  sec  r  in  ft 
s  s 

At  the  critical  radius  r^ 

icj(r*)  =  *** 19  ^  sec 

Using  Equations  (19)  and  (30)  the  ground  impulse  of  an  ideal  FAE 
in  air-MAPP  with  <£  =  0.  563  will  be 

I(r  )  =  9.  36  r  3  lb  sec  ;  r  <  r 

s  s  s  * 

T(r  =  3.06  r  3  -  11. 19  +  22.3 
s  s 

3 

=  (3.06r  +11.  lib  sec  ;  r  >  r* 

s  s 

The  variation  of  I(r  )  with  r  for  the  ideal  FAE  considered  here  is 

s  s 

shown  in  Figure  11. 

From  Equations  (33)  and  (35)  the  impulse  expressed  as  a  function 
of  time  becomes 

I(t)  =  C.  17  x  105  t)  lb  sec  ;  t  <  1.  88  x  10  ^  sec 

I(t)  =  (4. 60  x  10**  t3  +  19. 3)  lb  sec  :  t  >  1. 88  x  10  ^  sec 
The  variation  of  I(t)  with  t  is  shown  in  Figure  12. 


-  .  .  6.28(2.29)(10"3)(5.32)(103) 

ACJ  Fs  3.0 


0. 138  -  - 


2.  0(1. 30) (4. 69) 


39 


D.  DISCUSSION 


Simple  relations  are  presented  above  for  the  computation  of  the 
total  ground  impulse  generated  by  an  Idealized  FAE.  Given  the  detona¬ 
tion  properties  and  the  energy  released  by  the  secondary  blast,  the  rela¬ 
tions  presented  here  permit  computation  of  ground  impulse  with  minimal 
effort. 

The  idealized  FAE  involves  a  number  of  approximations.  The  self- 

(4  5) 

similar  strong  blast  wave  solution1  ’  '  is  used  to  represent  the  secondary 
blast.  The  perfect  gas  assumption  is  used  in  the  region  behind  the  lead¬ 
ing  shock,  back  pressure  is  neglected,  and  no  attempt  is  made  to  assess 

how  the  method  of  blast  initiation  affects  the  wave.  As  indicated  by 

(7) 

Brode,  Glass,  and  Oppenheim  ,  this  self-similar  blast  wave  solution 
provides  only  a  crude  approximation,  and  with  recent  advances  in  nu¬ 
merical  computation,  more  accurate  theoretical  solutions  have  become 
available.  However,  the  secondary  blast  energy  will  generally  be 
very  small  compared  to  the  total  energy  release  in  a  FAE,  i.e. ,  r*,- 
the  critical  radius  will  generally  be  much  smaller  Uian  the  fuel  cloud 
radius.  In  the  example  treated  above  r*  was  only  1.  54  feet.  Thus,  even 
though  the  self -similar  blast  solution  may  be  crude  it  has,  in  any  case, 
a  small  influence  upon  the  total  impulse  generated  by  the  FAE,  and  the 
simplicity  of  this  solution  makes  it  very  attractive  for  the  model  pre¬ 
sented  here.  It  is  encouraging  that  there  is  good  agreement  between  the 
idealized  FAE  and  the  tests  in  the  segmented  shock  tube  reported  below. 
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The  side  relief  due  to  the  inert  gas  bounding  the  upper  surface  of  the 
cloud  is  neglected  in  the  ideal  FAE.  The  expansion  wave  which  propagates 
from  the  inert  gas  into  the  combustion  products  behind  the  C-J  detonation 
will,  no  doubt,  result  in  an  appreciable  drop  in  the  ground  impulse  gener¬ 
ated  by  a  real  FAE .  The  influence  of  side  relief  is  currently  under 
investigation. 

Beyond  the  outer  radius  of  the  fuel  cloud  the  FAE  will  continue  to 

propagate  as  a  decaying  shock  wave  which  provides  a  further  contribution 

to  the  ground  impulse.  The  propagation  of  this  wave  has  been  examined 

(8) 

by  R.  Kiwan  by  replacing  the  C-J  detonation  by  a  spherical  piston  moving 

outward  with  constant  velocity  and  then  examining  the  nature  of  the  wave 

development  after  the  piston  stops.  The  numerical  solutions  obtained  by 

(8) 

Kiwan.  look  very  similar  to  the  self  similar  blast  wave  solution  and 
suggest  that  the  FAE  behavior  beyond  the  fuel  cloud  may  be  approxi¬ 
mated  by  a  blast  wave  with  suitably  chosen  initial  conditions.  This  aspect 
of  the  FAE  is  also  unde.1  study. 

The  secondary  blast  wave  is  generated  by  an  explosive  line  source 
in  this  ideal  FAE  and  so  has  cylindrical  symmetry.  In  an  actual  FAE 
it  is  more  likely  that  the  secondary  blast  will  originate  from  a  concen¬ 
trated  or  point  source  of  explosive  and  so  will  be  spherical.  The  details 
of  the  transition  from  spherical  to  cylindrical  symmetry  during  the  initia¬ 
tion  process  will  be  quite  complex  and, of  course, are  not  included  in  the 
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idealized  FAE.  Finally,  ground  shocks  and  reflection  of  the  FAE  shock 
*ront  from  the  ground  in  those  cases  when  the  fuel  cloud  is  not  immediately 
adjacent  to  the  ground  plane  are  neglected. 

The  sample  computations  and  the  detonation  data  presented  here  have 
been  determined  for  fuel  in  gaseous  form.  The  C-J  conditions  will  be 
modified  slightly  when  the  fuel  is  in  liquid  form.  This  problem,  which 
is  discussed  in  Reference  1,  is  under  investigation. 
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SECTION  III 


EXPERIMENTAL  RESEARCH 

A.  INTRODUCTION 

The  current  objectives  of  this  phase  of  the  research  have  been  to 
resolve  in  a  controlled  manner  overall  propagation  details  of  blast  initi¬ 
ated  unconfined  gas-phase  and  two -phase  detonations.  The  functional 
dependency  of  blast  wave  propagation  details  upon  time  and  radius  and 
of  detonation  wave  details  upon  time,  radius,  critical  threshold  energy, 
equivalence  ratio,  and  fuel  properties  has  been  sought.  To  achieve  these 
objectives  the  basic  experimental  hardware  previously  described ^  has 
been  utilized  with  some  additional  alterations. 

B.  EXPERIMENTAL  FACULTY 
1.  Drop  Production  Technique 

The  basic  method  for  generating  a  cloud  of  fuel  drops,  employed  from 
the  beginning,  has  continued  to  perform  satisfactorily.  However,  problems 
stemming  from  continuous  blockage  of  needles  by  combustion  products 
has  dictated  that  an  alternate  sub -manifold  design  be  chosen.  Shown  in 
Figure  13  is  one  such  assembly  of  the  new  design.  The  essential  improve¬ 
ment  of  this  design  lies  with  its  use  of  hypodermic  needles  which  may  be 
individually  replaced  in  the  event  of  blockage.  This  will  not  completely 
eliminate  the  needle  plugging  problem  but  will  provide  a  much  easier 
technique  for  correcting  for  it. 
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2-  Instrumentation 


Work  has  continued  with  improving  and  expanding  chamber  instrumen¬ 
tation.  Through  these  efforts  notable  improvements  in  data  regularity 
have  been  observed.  The  pressure  transducers  initially  employed  are  con¬ 
tinuing  to  give  satisfactoiy  data.  The  pressure  switches  described  in  the 
previous  report,  which  were  employed  as  time-of-arrival  devices,  have, 
however,  been  replaced  with  others  of  an  alternate  design.  Figure  14  is 
a  drawing  of  the  switch  now  being  used.  The  primary  advantages  of  this 
alternate  design  are  its  much  higher  reliability  and  ease  of  refurbishment. 
The  design  and  development  ol  heat  transfer  gages  were  initiated  during  this 
period.  Some  hardware  is  currently  being  tested,  but  a  usable  design  has 
not  yet  been  established. 

Likewise,  during  the  period  of  this  report  the  design  and  construc¬ 
tion  of  hardware  to  allow  an  optical  study  of  the  wave  processes  was  begun. 
Initially,  a  design  was  attempted  wlich  made  use  of  chamber  windows  of 
the  size  of  Hie  entire  chamber  side  plates.  This  design  met  with  little 
success.  Thus,  a  new  design  was  arrived  at  and  hardware  was  subsequently 
obtained. 

Figure  15  is  an  exploded  view  of  tho  new  chamber  side  plates  and 
Pyrex  windows  which  comprise  the  proposed  design.  The  practical  feasi¬ 
bility  of  this  hardware  shall  be  examined  lhrough  subsequent  tests. 
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Figure  14.  Current  Pressure  Switch  Design 


4  8 


3. 


Gas -Phase  Detonation  Facility 


Gas-phase  detonation  studies  were  carried  out  during  the  period 
encompassed  by  lh’s  report.  In  recent  years  some  investigations  have 
been  made  in  an  attempt  to  establish  the  dependency  of  r.Tical  threshold 
energy  for  detonation  initiation  as  a  function  of  the  concentration  of  MAPP 
(hydrocarbon  gas)-air  mixtures.  Among  these  investigations  were  he 
tests  performed  in  the  confinement  of  the  Crawshaw -Jones  Apparatu!^ 
the  bag  tests  of  Bei.edick  et  al*0^  and  the  bag  tests  of  Collins^1  ^  The 
investigations  performed  in  the  Crawshaw -Jones  Apparatus  reported 
functional  dependencies  of  detonation  limits  upon  initiation  energy 
which  largely  went  unsupported  by  the  latter  two  investigations.  The 
two  bag  test  studies  were  of  the  same  unconfined  nature  and  produced 
nearly  identical  results. 

It  is  apparent  that  the  degree  of  confinement  of  the  experimental 
apparatus  has  an  effect  upon  the  determined  detonation  initiation  limits. 
Hence,  it  was  of  much  interest  to  undertake  MAPP-air  detonation  tests 
in  order  to  establish  conclusively  the  unconfining  nature  of  the  sectored 
chamber  and  hence  its  utility  to  the  determination  of  fundamental  in¬ 
formation.  This  information  is  of  general  interest  but  also  of  particular 
interest  to  the  Air  Force  and  FAS  problem. 
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The  gas  used  in  these  tests  was  MAPP,  a  hydrocarbon  gas  com¬ 
prised  of  methyl  acetylene,  propane,  and  propadiene  gases.  Some 
results  of  this  study  are  presented  and  discussed  in  sub-Section  C 
of  Section  III. 

Some  modifications  to  the  basic  sectored  chamber  were  necessary, 
and  additional  support  hardware  was  required.  A  schematic  of  the 
gas -phase  detonation  apparatus  is  given  in  Figure  16.  Shown  in 
Figure  17  are  the  major  components  of  the  support  hardware.  The 
system  employed  used  a  MAPP-ait  reservoir  wherein  a  given  concentra¬ 
tion  of  MAPP  and  air  co».ld  be  pre-mixed.  The  resert'oir,  which  had 
a  volume  of  one  cubic  foot,  was  charged  to  60  psig  and  contained 
enough  mixture  to  allow  completion  of  up  to  fifteen  experimental  rims. 

This  provided  enough  MAPP-air  mixture  to  allow  completion  of  runs  at 
each  concentration  with  buv  a  single  charging  of  the  reservoir. 

Several  constraints  on  the  handling  of  MAPP  in  such  reservoir 
system  had  to  be  satisfied  prior  to  successful  operation.  One  such 
constraint  required  that  if  the  condensed  phase  was  to  be  avoided  the 
storage  pressure  of  MAPP  must  be  below  80  to  100  psig  for  a  temperature 
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Solenoid  Valve 


Gas  Detonation  Apparatus 


between  60  to  80°  F.  Hence,  the  reservoir  was  charged  to  60  psig  to  ensure 
the  MAPP  gas  and  air  would  mix  uniformly  <md  remain  so  throughout  the 
duration  of  the  runs.  An  additional  constraint  required  that  the  MAPP 
gas  cylinder,  containing  the  fuel  in  its  condensed  phase,  be  positioned 
in  its  inverted  position. 

When  vapor  is  removed  from  a  vessel  containing  a  liquid  hydrocarbon 

mixture,  a  simple  batch  distillation  without  rectification  ^called  weathering) 

occurs.  The  compositions  of  liquid  in  equilibrium  with  vapor  for  the 

weathering  of  a  stabilized  MAPP  mixture  have  been  determined  by  other 

/jl2  2.3  14) 

researcher^  *  ’  :  Their  results  demonstrate  the  potential  hazards 

with  the  weathering  of  MAPP  mixtures.  Hence  in  order  to  avoid  non- 
uniform  MAPP  compositions  and  unnecessary  safety  hazards  stemming 
ircm  unstabilized  MAPP,  liquid  MAPP  was  metered  from  its  inverted 
storage  cylinder.  Only  enough  liquid  was  metered  out  and  allowed  to  ex¬ 
pand  to  a  vapor  as  was  necessary  to  achieve  a  desired  partial  pressure. 
Many  precautions  were  taken  to  insure  that  the  desired  composition  was  ihal 
which  was  obtained  in  the  reservoir.  Samples  of  reservoir  chargings  ere 
taken  and  analyzed  by  gas  chromai ographic  techniques.  The  degree  c-f  sucee 
achieved  in  obtaining  desired  MAPP- air  concentrations  is  shown  in  Table  VI 

The  sectored  chamber  required  plugs  for  the  holes  in  the  top  plate, 
through  which  the  capillary  needles  for  the  two-phase  studies  normally 
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TABLE  VI.  MAPP  -AIR  CONCENTRATION0.  BY  V  OLUME  AT 

whic:t  tests  were  performed 


Desired 

Obtained 

Percent 

Difference 

9.  5 

9.7 

-2.1 

8.5 

8.5 

0.0 

7.0 

6.8 

2.8 

6.0 

6.0 

0.0 

4.5 

4.3 

4.4 

8.75 

3.8 

-1.3 

3.3 

3.3 

0.0 

3.0 

2.9 

_ 

3.3 
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penetrate.  The  chamber  also  required  a  fixture  at  the  downstream  end  to 
hold  a  mylar  diaphragm  which  is  replaced  for  each  run.  These  changes 
along  with  clay  packed  around  the  blasting  cap  leads  and  silicone  sealant 
in  strategic  locations  proved  to  lx1  satisfactory  enough  to  pull  the  chamber 
vacuum  down  to  within  0.  2  to  0.  5  in.  of  a  hard  vacuum. 

The  procedure  for  a  given  gas  detonation  run  was  as  follows: 

1.  Clean  chamber  of  combustion  products. 

2.  Insert  blasting  cap - Detasneet  initiation  charge,  secure  breech, 

and  clay  packing. 

3.  Secure  mylar  diaphragm. 

4.  Evacuate  chamber,  20  to  30  minutes. 

5.  Recharge  to  1  atm  with  current  MAPP-air  mixture. 

6.  Evacuate  all  plumbing  of  MAPP-air  mixture  from  chamber  to 
reservoir. 

7.  Initiate  timing  sequence  to  detonate  mixture. 

C.  RESEARCH  RESULTS  AND  DISCUSSION 

1.  Introduction 

Experimental  FAE  research  efforts  have  remained  focused  upon 
evaluation  of  the  two  limiting  cases  for  an  FAE  model  —  the  strong  cylin¬ 
drical  blast  wave  and  the  cylindrical  Chapman -Jouguet  detonation  wave- 
Table  VII  gives  a  summary  of  experimental  tests  performed  during  the 
period  of  this  report.  As  shown,  further  blast  wave  runs  were  made  with 
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TABLE  VII.  SUMMARY  OF  EXPERIMENTAL  TESTS  PERFORMED 


Run  No. 

Type 

Explosive 

Charge 

(grams) 

Comments 

100-135 

Blast  Wave 

1.0-2.  25 

Open  Breech 

200-211 

Detonation 

0. 5-3.0 

Kerosene  2-Air 

212-217 

Blast  Wave 

0.5 

Open  Breech 

218-282 

Detonation 

0.  5-3.  5 

Kerosene  2-Air 

283-286 

Blast  Wave 

1.0-1.75 

Open  Breech 

287-289 

Detonation 

1.  25 

Kerosene  2-Air 

290-309 

Blast  Wave 

1. 0-3.0 

Closed  Breech 

310-362 

Blast  Wave 

Small 

Detonators  a 

Open  Breech 

363-418 

Elast  Wave 

Small 

Detonators 

Closed  Breech 

419-559 

Detonation 

1 _ 

Small 

Detonator  -3.  5 

MAPP-Air 

aDuPont  electric  detonators  E-94  and  E-101  series  having  energy 
releases  smaller  than  generally  employed — Atlas  detonator. 


both  an  open  and  closed  breech  configuration  the  distinction  being  the 
presence  of  a  clay  plug  surrounding  the  blasting  cap  ignition  wires  as  they 
exit  from  the  breech  for  the  closed  breech  configuration.  The  detonation 
runs  were  composed  of  two  basic  types,  two-phase  Kerosene  2-air  and 
gas-phase  MAPP-air.  The  designation  Kerosene  2  is  intended  to  distinguish 
it  from  the  previously  employed  lesser  refined  fraction.  Kerosene  1. 

An  in-depth  analysis  of  the  •.<.  \st  wave  data  yielded  information  which 
helps  to  establish  the  nature  of  che  sectored  chamber  as  an  experimental 
model  of  cylindricallv  propagating  waves.  Two  basic  developments  followed 
from  this  analysis.  An  energy  efficiency  of  the  sectored  chamber  was 
derived  which  relates  the  effective  energy  release  with  the  calculated 
or  maximum  allowable  energy  released  in  the  chamber.  The  second 
development  stems  from  a  quantitative  comparison  of  experimental  blast 
wave  regression  models  with  strong  blast  wave  theory,  thereby  yielding 
a  behaviorial  estimate  of  wave  processes  in  the  sectored  chamber.  A 
presentation  and  discussion  of  two -phase  and  gas -phase  detonation  results 
follow. 

2.  Blast  Wave  Results 

Experiments  conducted  in  the  sectored  chamber  without  fuel  present 
have  continued.  In  order  to  make  meaningful  use  of  the  blast  wave  data  received 
from  these  tests,  it  was  necessary  to  arrive  at  a  standard  reduction  technique. 
Once  established,  this  technique  should  yield  self-consistent  results  useful 
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in  describing  the  behavior  of  wave  processes  in  the  sectored  chamber. 

The  reduction  technique  finally  chosen  begins  with  a  translation  of  the 
rough  radius-time  data  from  the  raster  scope  photographs  for  a  given 
energy  level.  This  rough  data  is  then  smoothed  using  a  least  squares 
polynomial  regression  to  the  model 

T  =  a01  +  a02  R  +  a03  R  ^ 

Examples  of  this  polynomial  smoothing  technique  are  given  in  Figure  18. 
Shown  here  are  radius-time  plots  for  three  energy  levels  1.0  gram,  1.75 
grams, and  3.0  grams  of  Detasheet.  Displayed  in  the  plots  are  the  actual 
rough  data  points  and  the  subsequent  least  squares  fit  of  Equation  (27)  to  the 
data. 

Equation  (27)  with  the  corresponding  calculated  coefficients  is  taken 
to  represent  the  experimental  data.  The  equation  is  then  used  to  compute 
values  of  t^as  previously  described  in  Reference  1.  Recall  that  t  is  the  value 
of  time  necessary  to  relate  the  experimental  data  to  the  imaginary  origin. 

All  time  values  of  the  specified  energy  level  are  then  adjusted  by  the 
computed  t  .  The  adjusted  data  is  then  subsequently  regressed  by  means 
of  the  least  squares  fit  to  the  following  experimental  models: 

T  =  a:  R  +  a2  R2 

T  =  b1  R2 
c2 

T  =  Cj  R 
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08) 

(29) 

(30) 


The  reason  for  selection  of  these  particular  models  ,'s  that  the  experi¬ 
mental  data  should  reflect  second  order  behavior  as  predicted  by  strong 
cylindrical  blast  wave  theory.  It  was  further  desired  to  verify  that  the 
comparison  be  suostantially  independent  of  the  choice  of  mathematical 
regression  models.  Since  the  data  has  already  been  adjusted  to  the  imag¬ 
inary  origin,  no  leading  constants  are  necessary.  A  numerical  check 
was  performed  to  examine  just  how  clcse  to  the  origin  the  adjusted  lata 
curves  come.  For  this  purpose  the  adjusted  data  was  again  fit  to  the 
model  of  Equation  (30)  and  new  tQ's  computed.  This  check  in  all  cases 
suggested  that  an  iterative  process  seeking  t^  was  not  necessary  as  the  t 
in  the  second  run  was  within  a  small  number  of  being  zero.  Tnis  entire 
data  reduction  technique  was  repeated  for  each  energy  level  as  characterized 
by  a  given  detonator -Detasheet  combination. 

During  the  course  of  these  tests  and  following  data  analysis  it  was 
found  that  the  energy  release  experienced  in  the  sectored  chamber  was 
less  than  that  which  may  be  predicted  from  a  known  quantity  of  Detasheet. 
This  fact  is  quite  justifiable  on  the  basis  that  non-idealities  exist  in  the 
physical  model  of  a  cylindrical  blast  or  detonation  wave.  The  significant 
deviations  from  an  ideal  blast  wave  model  are: 

(1)  The  existence  of  a  physically  limiting  chamber  suggests  losses 
to  the  walls. 

(2)  Initiation  energy  release  is, in  fact.not  instantaneous  as  assumed 
by  theory.  It  was  estimated  to  be  of  the  order  of  10  psec. 

(3)  Initiation  energy  is  not  released  uniformly  along  a  line  source - 
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The  numerical  values  of  energy  efficiencies  were  arrived  at  through  use 
of  experimental  data  regressed  by  the  model  of  Equation  (29)  and  the 
general  radius-time  form  of  strong  cylindrical  blast  wave  theory 
given  in  Equation  (31) 


_  -i  /  9  9 

t  =  (Eo/Pl)  *'  V 


or 


(31) 


t  =  dj  r 


Since  Equations  (29)  and  (31)  are  of  the  same  general  form,  any  dif¬ 
ferences  in  radius -time  behavior  between  the  two  for  a  given  process 
relates  directly  to  deviations  in  their  constant  "coefficients  b^  and  dj. 
It  is  clear  then  that  b^  reflects  actual  behavior  while  d^  reflects 
theoretical  behavior  in  the  sectored  chamber.  This  fact  can  be  stated 
by  the  representation  of  b^  and  d^  as 


bi  = 


dl  = 


(E/Pi)'172 


Jregression 


(E0/Pj)"1/2 


calculated 


Then  E  and  E  are  given  bv 
o  o  , 

reg  calc 


E0  =P1/b1 
reg 


En  =  p,/d2 

0  ,  A  1 

calc 


(32) 


(33) 
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An  efficiency  factor  can  now  be  constructed 


r?  =E  /E 
e  c  o  . 
reg  calc 


(34) 


where  the  energies  are  defined  bv  Equations  (33).  Thus  En  follows 

J  ureg 

directly  from  a  mathematical  regression  to  experimental  data,  while  E. 


calc 


is  derived  from  knowledge  of  th2  detonator  -Detasheet  combination.  For 
the  combination  of  an  Atlas  blasting  cap  yielding  1110.  0  ft-Ibf  energy  and 
X  grams  of  detasheei  yielding  2911.  0  ft-lbf/gram,  E  in  the  sectored 


calc 


chamber  would  become 


Eq  =  (12.  0/2.  05)(360/ 20) (1 1 10.  0  +  2911  X)  ft-lbf/ft  (35) 

calc 

Energy  efficiencies  can  then  be  computed  for  a  range  of  energies  for 
which  experimental  blast  wave  data  was  taken.  Figure  19  gives  the  results 
of  this  computation  for  both  open  and  closed  breech  configurations.  It  was 
deemed  necessary  to  assess  the  degree  which  closed  breech  wave  behavior 
deviated  from  the  open  breech  case.  The  former  configuration  was 
employed  in  the  MAPP-air  gas-phase  detonation  tests.  The  curve  dis¬ 
played  is  a  second  order  polynomial  regression  to  all  the  data  and  was 
used  to  describe  the  data  in  subsequent  computations.  The  closed  breech 
data  fell  within  the  data  scatter  and  hence  no  definite  trend  of  its  own 
was  determined. 

In  the  interest  of  examining  how  closely  cylindrical  behavior  was 
modeled  in  the  sectored  chamber,  experimental  data  was  compared  with 
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Figure*  19.  Sectored  Chamber  Energy  Efficiency  as  a  Function  of  Calculated  Energy  Release. 


strong  blast  wave  theory.  A  quantitative  comparison  between  experiment 
and  theory  was  achieved  by  making  use  of  the  experimental  regression 
models  of  Equations (28\  (29)and  &Dan G  the  following  expressions  for 
cylindrical  blast  theory 


1/  4 

-O' 


1/2 


(36) 


■wUO 


1/4 


-1/2 


(37) 


The  energy  used  in  the  above  theoretical  expressions  was  that  given 
by  Equation  (34)and  repeated  here  as 


E  J  ,  =  E  .xi?  (ft-lbf/ft) 
cylindrical  calc  e 


(38) 


Figures  20,  21, and  22  reveal  the  radius-time  and  Mach  number-radius 
behavior  for  0. 0,  1.  5, and  3.  0  grams  of  Detasheet, respectively.  Shown 
in  these  figures  are  the  three  experimental  regression  models  as  com¬ 


pared  with  strong  cylindrical  blast  wave  theory. 

All  the  curves  of  Figures  20,  21, and  22  are  self-consistent  and 
substantially  in  agreement  about  major  behavioral  trends.  However , 
Equation  (29),  the  form  for  ideal  strong  cylindrical  blast  wave  behavior, 


is  seen  to  give  the  much  better  agreement  between  experiment  and  theory. 
In  part  this  is  to  be  expected  because  the  effective  energy  release  was 
obtained  by  assuming  such  behavior.  On  the  other  hand,  the  adjusted 
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Chamber  Radius,  R  (in.) 


Figure  20(b).  Mach  Number  versus  Radius  Behavior  of  Strong  Blast 
Wave  Theory  and  Experimental  Regression  Models  for  0. 0  gram  o. 

Datasheet 


Imcmtul  Regression  Models  for  1. 5>  grams  of  Dotnshool 


e  22(a),  Radius  versus  Time  Behavior  cf  Strong  Blast  Wave  Theory  and 
Experimental  Regression  Models  for  3.0  grams  of  Detasheet  fC\ 


Chamber  Radius,  R  (in.) 


Figure  22(b).  Mach  Number  versus  Radius  Behavior  of  Strong  Blast 
Wave  Theory  and  Experimental  Regression  Models  for  3.0  grams  of 

Detasheet  fC'. 
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time  coordinate  was  arrived  at  by  a  different  technique. 

Inasmuch  as  the  theoretical  curve  involves  use  of  an  experimentally 

determined  efficiency  factor,  the  influence  of  this  value  on  predictions 

was  examined.  Various  polynomial  fits  to  the  efficiency  data  of  Figure  19 

were  assumed  and  showed  little  influence  on  the  value  of  77  for  a  given 

energy  level.  Further,  variations  in  f]  were  assumed  and  the  changes 

e 

in  theoretical  blast  wave  arrival  times,  wave  velocity,  and  wave  Mach 
number  were  noted.  A  10  percent  variation  in  r)Q  led  to  a  3  to  4  percent 
variation  in  the  latter  quantities. 

It  is  concluded  that  the  experiments  do  exhibit  the  desired  cylindrical 
behavior  and  that  this  characteristic  is  taken  on  rather  early,  say  by  a 
radius  of  6  inches. 

It  is  anticipated  that  the  blast  wave  phase  of  this  research  will  not 
be  terminated  completely.  Periodically  additional  runs  will  be  made  to 
gain  a  good  statistical  sample  for  establishing  the  energy  efficiency. 

Mso  a  comparison  of  blast  wave  pressures  and  possibly  heat  transfer 
as  functions  of  radius  is  planned. 

3.  Two-Phase  Detonation  Results 

Experiments  have  been  conducted  in  monodisperse  sprays  of  a  highly 
refined  fraction  of  Kerosene,  subsequently  referred  to  as  Kerosene  2, 
the  oxidizer  was  air  at  atmospheric  conditions.  At  this  time  the  controls 
exercised  upon  the  conditions  of  the  detonation  runs  have  been  limited  to 
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systematic  variations  in  initiation  energy,  Eq,  at  a  fixed  global  equivalence 
ratio.  ,  for  a  given  fuel.  Variations  in  to  investigate  the  influence 
of  this  parameter  upon  wave  properties  are  easily  accommodated  and 
this  is  planned  for  future  experiments.  Presently  detonation  runs  have 
been  made  at  $  ~  0.  63.  This  figure  was  calculated  on  the  basis  of  an 
experimentally  determined  mean  drop  size  of  380  microns.  A  5  percent 
variation  in  mean  drop  size  by  these  calculations  produces  a  10  percent 
variation  in  <£_. 

The  liquid  fuels  used  to  date  and  those  with  which  further  experiments 
are  planned  are  given  in  Table  Vin,  along  with  their  more  common  physical 
properties.  As  expected,  the  properties  of  the  two  kerosene  fractions 
are  very  much  alike,  as  are  the  properties  of  the  two  nitropropane  types. 

In  view  of  this,  1 -Nitropropane  may  be  dropped  from  the  experimental 
schedule.  The  Gordon-McBride  NASA  program  which  has  proved  so 
useful  in  the  past  has  been  used  to  derive  equivalent  gas  phase  detonation 
properties.  Plotted  in  Figures  23,  24,  and  25  are  these  properties  for 
Kerosene  2,  2 -Nitropropane  and  Propyl  Nitrate, respectively .  A  cursory 
comparison  of  the  curves  reveals  behavioral  trends  which  should  be  inter¬ 
esting  ne  experimentally.  Detonation  velocities  for  the  two  latter 

fuels  mentioned  above  tend  to  effectively  plateau  for  a  large  range  of  equiva¬ 
lence  ratios.  Additionally,  the  behavior  in  the  pressure  ratios  of  these 
same  two  fuels  follows  trends  considerably  different  from  the  fuels  which 
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Figure  25(b).  Propyl  Nitrate -Air  Gao  Detonation  Parameter. 


TABLE  VIII.  COMPARISON  OF  PHYSICAL  PROPERTIES  OF  HYDROCARBON  FUELS 


Equivalence  Ratio 
(O0  basis) 


are  currently  being  examined.  Pressure  ratio  behavior  of  this  type  is 
partially  due  to  the  increased  amounts  of  oxygen  present  in  the  reaction. 

Reduction  of  the  experimental  detonation  data  proved  it  to  be  very 
self-consistent  and  in  substantiating  agreement  with  that  reported  upon 
previously.  Figures  26  and  27  are  characteristic  radius-time  plots  of 
detonation  data  at  varying  energy  levels  for  Kerosene  1  and  Kerosene  2. 
These  plo*s  demonstrate  the  general  form  of  the  rough  data  as  well  as  the 
currently  employed  fourth  order  polynomial  fit  to  this  data.  The  reason 
for  selecting  a  fourth  order  polynomial  was  simply  that  it  tended  to  reflect 
wave  properties  sufficiently  well: 

1.  Slope  of  r-t  plot  beyond  rt  is  generally  constant  for  appropriate 
energies,  if  low  frequency  polynomial  generated  oscillations 
are  ignored. 

2.  Constant  detonation  velocity  is  predicted  for  energies  above 
the  critical  threshold  energy. 

3.  Decaying  detonation  velocity  is  predicted  for  energies  below 
the  critical  threshold  energy. 

4.  Critical  radii  is  predicted  with  remarkable  accuracy. 

5.  Detonation  velocity  similarly  predicted  with  remarkable  ac-uracy. 
This  technique,  however,  does  possess  many  disadvantages,  none  the 

least  of  which  is  a  lack  of  sensitivity  to  variations  in  wave  velocity.  Con¬ 
sequently,  work  on  an  alternate  detonation  data  regression  mode1  was 
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Figure  27(g)*  Kerosene  2-Air  Detonation  Data,  2.5  grams 


Tlmt/100  n9ec 

Figure  27(h).  Kerosene  2-Air  Detonation  Data,  3.0  grams 


Tlme/100  tine c 


Figure  27(i).  Kerosene  2-Air  Detonation  Data,  3.  5  grams 


undertaken.  This  regression  model,  which  is  not  yet  completed,  makes 


use  of  the  method  of  LaGrange  multipliers  to  minimize  the  standard 
deviation  of  a  least  squares  curve  fit  of  rough  experimental  radius-time 
data  to  the  following: 

9, 

for  R.  <  r* 


T.  =  a  +  a1  R.  +  Q'nR.  +  €■ 

1  0  1  1  2  l  1 


and 


(39) 


T.  =  (3  +  P*  R.  +  e.  for  R.  >  r^ 

l  Ho  1  1  i  l  —  * 

The  problem  is  then  to  compute  the  constant  coefficients  aQ,  o-p  a ^  PQ. 

|3p  and  the  critical  radii  r*  in  order  to  minimize  the  standard  deviation 

subject  to  the  constraints  that  at  R^  =  r* 


and 


a q  +  +  <*2T*  ~  ^ 1 


+  2^2  =  3-^ 


(40) 


The  primary  advantages  to  this  model  are  the  greater  reliability  in  the 
derived  critical  radii  values  and  the  well  behaved  continuous  nature  of 
the  final  function  T  =  T(R).  The  latter  point  will  allow  for  a  more  repre¬ 
sentative  continuous  record  of  wave  propagation  velocity. 

Figure  28  demonstrates  the  current  technique  for  obtaining  detonation 
wave  propagation  velocity  and  critical  radii.  Displayed  in  Figure  29  and 
Table  IX  are  the  basic  results  for  Kerosene  1.  Figure  29  is  a  plot  of 
two-phase  detonation  velocity  as  a  function  of  equivalence  ratio  and  drop 
size.  The  predicted  theoretical  two-phase  velocity  was  determined  by 
modifying  the  equivalent  gas  phase  velocity  by  a  mean  velocity  deficit. 
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Figure  28(a).  Data  Reduction  Technique,  1.  5  grams 


Figure  28(b).  Data  Reduction  Technique,  2.  5  grams 
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29  .  Spray  Detonation  Velocity  for  Kerosene  1  as  a 
Function  of  Equivalence  Ratio  and  Drop  Size. 


TABLE  IX.  EXPERIMENTAL  TWO-PHASF  DETONATION  RFSULTS  FOR  KEROSENE  \. 


Recall  that  velocity  deficit  is  a  function  of  drop  size,  chamoer  hydraulic 
radius,  fuel,  and  oxidizing  medium.  The  mean  or  average  velocity  deficits 
used  in  the  determination  of  kerosene  spray  velocity  were  computed  to  be 
1.44%.  2.  88%  and  4.  32%  for  200  p,  400pand  600p  drop  sizes, respectively. 
The  experimental  values  of  spray  velocity  and  critical  radii  are  given.  The 
results  appear  quite  satisfactory.  Similarly,  Figure  30  and  Table  X  give 
the  detonation  results  for  Kerosene  2.  The  comparison  between  experi¬ 
ment  and  theory  proved  to  be  even  more  satisfactory  than  with  Kerosene  1. 
Definite  trends  in  this  comparison  are  now  beginning  to  show  up.  Basically, 
the  difference  between  experiment  and  theory  becomes  wider  as  the  explo¬ 
sive  charge  is  decreased  toward  the  ignition  threshold  energy.  This  holds 
true  for  velocity  as  well  as  critical  radii  data  for  both  kerosene  fractions. 

Figure  31  displays  in  a  more  transparent  fashion  the  fundamental 
wave  behavior  of  a  Kerosene  2-  air  detonation  as  a  function  of  chamber 
radial  distance  and  explosive  charge.  The  continuous  and  expected  decay 
patterns  of  the  curves  up  to  the  critical  radii  are  clearly  noticeable.  The 

trends  depicted  here  are  in  contrast,  but  not  necessarily  in  disagreement, 

(l) 

with  those  reported  upon  in  Figure  48  of  the  previous  annual  i  eport  . 

The  dissimilar  decay  trends  reported  upon  therein  suffered  by  virtue  of 
the  fact  that  no  mathematical  regression  model  was  available  with  which  to 
describe  the  experimental  detonation  data.  Consequently,  velocities  were 
obtained  by  taking  slopes  of  radius-time  data  curves  by  observation.  Such 
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TABLE  X.  EXPERIMENTAL  TWO-PIIASE  DETONATION  RESULTS  FOR  KEROSENE  2. 
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energy.  During  the  construction  of  these  curves  it  was  observed  that 
the  experimental  blast  wave  rate  of  decay  quickened  with  increasing 
energy.  Hence,  the  deviation  of  the  detonation  from  that  of  strong 
blast  wave  M-r  behavior  was  increasingly  more  evident  as  energy  was 
increased.  This  deviation  is  primarily  due  to  the  fact  that  while  blast 
wave  velocity  decays  faster  with  increasing  energy,  the  critical  radii 
for  detonation  increases.  It  was  further  observed  that  at  the  critical 
threshold  energy  for  Kerosene  2  (approximately  1.  5  grams)  the  experi¬ 
mental  blast  wave  datum  curve  nearly  passed  through  the  critical  radii- 
detonation  Mach  number  intersection.  It  would  appear  from  this  crude 
representation  that  the  larger  the  blast  energy  used  to  detonate  a  fuel, 
the  quicker  the  attendent  blast  decay  and  the  sooner  the  energy  released 
by  the  detonation  process  has  a  modifying  effect  upon  this  decay.  All  this 
occurs  subject  to  the  additional  requirement  that  the  proper  critical  radii- 
detonation  Mach  number  constraint  is  satisfied  at  each  energy  level. 

The  general  form  of  the  experimental  M-r  curves  of  Figure  31  then 
presents  the  picture  of  a  wave  process  initially  dominated,  but  slightly 
modified,  by  the  generally  characteristi  decay  of  a  strong  cylindrical 
blast  to  an  energy  dependent  critical  distance.  At  *his  critical  point  an 
abrupt  transition  io  a  detonation  dominated  wave  process  occurs.  A 
dependency  of  the  final  detonation  velocity  upon  the  initiation  energy 
appears  to  show  up  to  some  degree  in  these  curves.  A*  this  time,  however, 
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such  a  dependency  can  only  be  conjecture  in  view  of  the  error  tolerances 
associated  with  the  data  reduction  technique  used  to  obtain  these  curves. 
Many  such  dependencies  should  be  more  thoroughly  established  with  the 
subsequent  use  of  the  detonation  regression  model. 

Plotted  in  Figure  32  is  critical  radii  as  a  function  of  calculated 
energy.  Two  analytic  curves  appear  with  experimental  points  for  both 
kerosene  fractions  superimposed.  The  upper  curve  is  critical  radii 
computed  as  a  function  of  energy  by 


where  100  percent  of  E0  calculated  is  used.  The  lower  curve  is  critical  radii 
computed  by  Equation  (41)  where  Eo  is  altered  by  blast  wave  energy  efficiency 
discussed  earlier.  The  results  are  most  striking,  since  a  remarkable 
agreement  occurs  between  an  analytically  determined  r*  and  experimental 
values.  The  same  value  of  P^/P^  =  12.  55  was  used  to  compute  Q  for 
both  kerosene  fractions,  thus  making  possible  the  representation  of 
experimental  r*  values  on  the  same  plot. 

4.  Gas-Phase  Detonation  Results 

Experiments  have  been  conducted  in  gaseous  MAPP-air  mixtures. 
Experimental  runs  were  made  in  these  mixtures  for  a  constant  pressure 
of  1  atmosphere  and  atmospheric  temperatures  ranging  from  19.  0°C  to 
25.7°C.  The  controls  exercised  upon  the  conditions  of  the  detonation 
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Critical  Ignition  Distance  as  a  Function  of  Energy  for  Kerosene  1  and  Kerosene  £ 


runs  were  over  a  broader  range  than  those  presently  exercised  upon  the 
two-phase  detonation  runs.  Experimental  data  on  wave  propagation  as 
a  function  of  radius  was  obtained  for  a  systematic  variation  of  initiation 
energy  at  each  of  eight  different  equivalence  ratios  as  given  in  Table  VI. 

It  was  deemed  necessary  tn  establish  the  composition  of  the  MAPP- 
gas  used  for  the  detonation  tests.  Consequently, a  gas  chromatographic 
analysis  was  performed  upon  several  samples  of  the  pure  gas  as  well  as 
upon  the  eight  MAPP-air  mixtures.  The  analysis  was  performed  upon  a 
Varian  90-C  gas  chromatograph  with  the  conditions  specified  in  Table  XI. 

Figures  33  and  34  give  examples  of  gas  chromatograms  produced 
by  the  analysis.  As  is  apparent  from  the  chromatograms,  separation  of 
the  constituents  was  very  well  defined  despite  the  fact  that  they  were  very 
nearly  of  the  seme  molecular  weight.  Figure  33  is  a  chromatogram  for 
pure  MAFP,  while  Figure  34  is  for  a  10  percent  MAPP-air  mixture.  The 
jump  in  the  base  line  between  air  and  propane  is  due  to  a  scale  change 
necessary  to  keep  the  traces  on  scale.  Table  XII  lists  the  significant 
properties  of  the  MAPP  gas  used  during  the  tests.  Once  established 
these  properties  were  used  as  input  data  to  obtain  characteristic  gas 
detonation  parameters  from  the  previously  mentioned  NASA  program. 
Given  in  Figure  35  are  plots  of  these  parameters.  An  anomalous  density 
ratio  behavior  shows  up  for  equivalence  ratios  nearing  3.  0.  The  primary 
reason  for  this  is  due  to  the  appearance  of  solid  carbon  in  the  combustion 
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Figure  33 .  MAPP  Gas  Chromatogram. 


Acetylene 


Figure  34.  MAPP -Air  Gas  Chromatogram  (10%  MAPP) 


MAPP-Air  Gas  Detonation  Properties. 


Figure  35(b).  MAPP -Air  Gas  Detonation  Properties. 


TABLE  XI.  GAS  CHROMATOGRAPH  ANALYSIS  CONDITIONS 


Column: 

Length 

10  ft 

Inside  diameter 

1/4  in. 

Composition 

QF-1  30%  on  60-80  chromosorb 

Settings: 

Column  temperature 

55°F 

Detector  temperature 

150°F 

Injector  temperature 

30°F 

Helium  flow 

11  ml/min 

Filament  current 

100  milliamperes 
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TABLE  XII.  MAPP  GAS  PROPERTIES 


Molecular  Weight 

41.0 

Density  (slugs/ ft**) 

0.1126 

Stoichiometric  equivalence 

ratio  (0^  basis) 

0.302 

Composition 

Methyl  Actylene 

51.0% 

Propane 

26.0% 

Propadiene 

23.0% 

Appearing  in  Figure  36  is  a  sample  series  of  radius-time  plots  of 
rough  detonation  data  for  concentration  of  9.  7  percent  MAPP  (by  volume) 
in  air.  The  threshold  ignition  energy  for  this  concentration  was  selected 
to  be  1.6  grams  Detasheet.  It  is  at  this  energy  where  the  curve  first 
begins  to  suggest  a  constant  slope,  implying  constant  propagation 
velocity  beyond  a  given  radii.  A  similar  series  of  curves  were  generated 
for  each  of  the  otner  seven  concentrations  to  establish  the  dependence  of 
critical  threshold  energy  upon  MAPP  concentration.  This  dependence  is 
displayed  in  Table  XIII,  and  it  is  plotted  in  Figure  37.  The  characteris¬ 
tically  narrow  U-shape  of  the  MAPP -air  detonation  curve  has  been 
satisfactorily  established.  The  limits  of  detonability  found  for  the  MAPP 
gas  used  were  from  2.  9  percent  to  10.  5  percent  by  volume.  The  rich 
limit  is  an  extrapolation  of  data  taken  up  to  MAPP  concentration  of  9. 7 
percent  by  volume  and  represents  a  reasonable  extension  of  the  data. 

It  is  of  interest  to  examine  these  detonation  limits  in  light  of  other 
recent  experimental  results.  Table  XIV  represents  a  comparison  of 
detonation  limits  on  a  percent  volume  basis  between  four  separate  experi¬ 
mental  studies.  By  using  the  Crawshaw -Jones  Apparatus,  it  was  found 
that  detonation  limits  of  MAPP-air  mixtures  widened  with  increasing 
initiator  energy.  Yet  the  bag  tests  of  Benedick  et  al.  and  Collins^1" 
did  not  confirm  this  functional  dependency.  Similarly,  the  present  study. 


Threshold  Energy  x  10 “J  (fl-lbf/ft) 


TABLE  X1U.  CRITICAL  ENERGY  THRESHOLD  FOR 
MAPP- AIR  MIXTURES  IN  THE  SECTORED  CHAMBER 


<3? 

Percent  MAPP  in 

Air  by  Weight 

Percent  MAPP  in 
Air  by  Volume 

Effective  Initiator 
Energy 

mmm 

2.  20 

13.33 

9.7 

3875.3 

408, 664 

1.90 

11.69 

8.5 

2243.0 

236, 190 

1.49 

9.43 

6.8 

1011.9 

106, 553 

i.30 

8. 36 

6.0 

793.4 

83,  545 

0.90 

5.92 

4.3 

647.9 

68,  234 

0.78 

5.12 

3.8 

1144.1 

120,471 

0.59 

4.40 

3.3 

2039. 2 

214,731 

0.56 

4.00 

2.9 

4160.7 

438,118 
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TABLE  XIV.  DETONATION  LIMITS  OF  MAPP-AIR  MIXTURES  BY  VOLUME 


1 — - — - 

Method 

Initiator 

Lower 

Limit 

Upper 

Limit 

Crawshaw-Jones  Apparatus 

1  gram,  PETN 

4.1 

7.8 

Crawshaw -Jones  Apparatus 

10  grams,  PETN 

2.4 

13.7 

Crawshaw-Jones  Apparatus 

100  grams,  PETN 

- 

>  30 

Bag  Test 

800  grams,  C-4 
(672  grams,  PETN 
equivalent) 

2.9 

10.2 

Bag  Test 

1 

386  grams,  PETN 

2.9 

9.1 

Sectored  Chamber 

a 

2  grams,  Detasheet 
'O’  (1.  57  grams, 
PETN  equivalent) 

2.9 

10.  5a 

Extrapolated  from  test  results  taken  up  to  9.  7 percent  by  volume. 
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using  initiator  energies  two  orders  of  magnitude  smaller  than  the  latter 
two  studies,  tends  to  cast  doubt  upon  such  a  dependency.  There  is  ob¬ 
viously  a  need  for  further  work  on  this  subject. 

A  further  comparison  was  performed  to  examine  the  critical  threshold 
energy  limits  of  the  present  study  relative  to  those  of  recent  AFATL  bag 
tests  by  P.  Collins^  \  Figure  38  is  a  plot  of  nondimensional  critic^ 
threshold  energy  as  a  function  of  MAPP  concentration  by  weight  for  the 
two  mentioned  studies.  The  nondimensional  energy,  E,  was  arrived 
at  by  dividing  all  energies  for  a  given  study  by  an  energy  selected  from 
this  data  as  the  standard.  The  standard  for  each  study  was  the  value 
of  critical  energy  corresponding  to  11.7  percent  MAPP  by  weight. 

This  standard  is  not  absolute,  but  rather  it  was  conveniently  chosen  to 
demonstrate  relative  trends  in  critical  energy  since  data  was  obtained 
in  both  studies  at  this  MAPP  concentration.  The  comparison  reveals  the 
present  study  produced  data  suggesting  a  slightly  broader  and  slightly 
shallower  characteristic  threshold  curve.  The  broader  detonability 
limits  in  the  present  study  are  supported  by  the  fact  that  the  MAPP  gas 
used  had  a  higher  percent  of  methyl  acetylene  present  as  compared  to 
the  MAPP  used  in  the  bag  tests.  The  composition  of  the  MAPP  used  in 
the  bag  tests  is  reported  to  have  approximately  the  composition  37  percent 
methyl  acetylene,  25  percent  prodadiene,  20  percent  propane,  9  percent 
C^-carbon  compounds  (mostly  n-butane)  by  volume. 


Ill 


Figure  38.  Comparison  of  Current  MAPP-Air  Detonation 
Initiation  Limits  with  AFATL  "Bag"  Test  Results. 
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It  is  apparent  from  the  foregoing  discussions  that  the  experimental 
phase  of  this  research  continues  to  lend  encouragement  to  the  prediction 
of  gas -phase  and  two-phase  detonation  wave  properties  and  threshold 
energy  levels  required  for  detonation  initiation. 
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For  the  blast  wave  solution  V  lies  in  the  range 


<  V  < 


(v  +  2)y  -  —  (v  +  2) (y  +  1) 


(I  -7) 


The  value  of  V  =  4/(v  +  2)(y  4  1)  corresponds  to  a  point  immediately 
behind  the  shock  wave  and  the  point  V  =  2/{v  +  2)y  corresponds  to  the 
singular  point  at  the  center  of  the  blast  wave. 

Equations  (I  -2)  to  (  1-6)  have  been  programmed  for  the  computer 
to  produce  the  variation  of  v-Vg.  p/p9«  p/pp,  and  T/Tg  as  a  function 


116 


of  r/r  .  The  calculation  is  started  at  V  -  4/(y  +  2 )(y  +  1)  and  V  is 

allowed  to  approach  the  blast  center  value  of  2,  (c  ^  2)") .  As  a  check 

(4) 

several  cases  computed  by  Sedov  were  also  determined  using  the  blast 
computer  program  and  exact  agreement  was  found. 

A  separate  computer  program  was  developed  to  compute  the  param¬ 
eters  p(y,  v)  from  Equation  (15)  and  ot(y,v  )  defined  by 


a(y,  v)  =  O' 


\  \  R  V2\V+l  dX  +  - 

2  J  y 


1 

hi 


pa"*1  dA 


(1-8) 


It  is  important  to  recognize  that  p/pg  and  the  similarity  variable  P(A) 
are  connected  by  the  relation 


P=/JL 


\2/{v  +  2 Hy  +  1) 


The  above  computer  programs  and  sample  outputs  are  available  on 
request.  Typical  variations  of  p/pg,  v/vg,  anc*  T/Tg  with  the  ratio  of 
specific  heats  y  and  with  the  geometry  parameter  v  are  shown  in  Figures 
1-1  to  1-6. 
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Dimensionless  Pressure  vs  Distance 


APPENDIX  n 


CALCULATION  OF  THE  GROUND  IMPULSE 
FROM  A  CBAPMAN-JOUGUET  DETONATION 


As  indicated  in  Reference  1  the  conservation  equations  describing  the 
self-similar  flow  behind  a  Chapman -Jouguet  Detonation  can  be  reduced  to 
the  single  ordinary  differential  equation 


dz  z  [2(V-1) 2  +  (u  -  l)(r2-l)V(V-l)  -  2z] 
dV  V[(V-1)2  -  V7-] 

where  V  is  a  dimensionless  velocity  defined  by 


(II-l) 


v 


( II  -  2) 


and 

z  =  r  i  ;  p  =  pAr2/t2)'P 

R  1  ( II-  3) 

P  =  PjR 

y  is  the  ratio  of  specific  heats  of  the  combustion  products.  Integration  of 
Equation  (Il-l),  which  must  be  carried  out  numerically,  is  the  key  problem 
in  determining  the  self-similar  flow  behind  a  C- J  detonation.  Once  the  re¬ 
lation  between  z  and  V  is  determined,  the  variation  of  V  and  z  with  the  di¬ 
mensionless  radius  A  =  (r/rj  can  be  found  by  integrating  the  equation 


d(ln  A) 
dV 


z  -  (v-D2 

V[(V-1)2-!/z] 


(II  -4) 


The  flow  behind  the  detonation  is  assumed  to  be  isentropic  so  that  z 
and  R  are  related  by 


(y2-« 


R 


(II— 5) 


lv 


Zg  and  R^  are  the  values  of  z  and  R  immediately  benind  the  detonation  front 
and  for  a  C-J  detonation  are  given  by 


y2  Z1  -1 

r„  =  [  +-^3 

2  L  y2+ 1  J 


y2  Z1  2 

z  _  - £ —  q  +-^ 

2  ,  ^2'  y/ 

(y„+l)  1 


z^,  the  value  of  z  immediately  ahead  of  the  detonation  is  given  by 


(n-6) 


M 


z  =  y 


'i^A2) 


1  Rh 


(1.0) 


(II-7) 


since  R.  =  1.0.  Since  r  =  Ct  for  a  C-J  detonation 
1  s 

ri  — 


zi  = 


m: 


D 


(11-8) 


The  behavior  of  the  C-J  solution  in  the  z-V  plane  as  determined  by 
Equation  ( II- 1)  has  been  discussed  in  Reference  1,  and  is  also  indicated  in 
Figure  ( IX— 1)  below-  Ahead  of  the  detonation  front  v  and  hence  V  =  0  while 


p  =  p.^  and  p  =  Pj.  Thus,  the  variable  z  becomes 
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Figure  II-l.  The  v-z  Plane  for  a  C-J  Detonation. 


The  undisturbed  fuel -oxidizer  mixture  at  r  — >  oc  thus  corresponds  to  the 

point  z  =  0,  V  =  0,  while  the  point  immediately  ahead  of  the  detonation  front 

_2 

corresponds  to  z  =  z^  =  ;  V  =  0.  There  is  then  a  discontinuous  jump 

across  the  detonation  front  to  the  point  ( z V^)  immediately  downstream  of 

2 

the  detonation.  For  a  C-J  detonation  (z^,  Vg)  lies  on  the  parabola  z  =  (1-V) 
the  locus  of  points  where  the  velocity  is  sonic  with  respect  to  the  detonative 
discontinuity. 

The  detonation  front  is  followed  by  an  isentropic  expansion  through 
which  the  velocity  of  the  combustion  products  drops  to  zero  at  the  boundary 
of  a  stationary  core  region,  which  corresponds  to  the  point  A:  z  =  1,  V  =  0. 

The  point  A  is  a  singular  point  which  can  be  shown  to  be  a  node, and  it  is  the 
transition  from  ( z 2, 
tion  of  Equation  (H-i).  In  the  physical  plane  the  point  A  moves  radially 
outward  with  the  speed  of  sound  and  corresponds  to  the  characteristic  sepa¬ 
rating  the  stationary  core  from  the  expansion  behind  the  detonation.  The 
stationary  core  is  represented  by  the  line  V  =  0  extending  from  z  =  1  to 
z  =  cc  which  corresponds  to  the  detonation  center  r  =  0. 

At  the  singularity  A,  Equation  (n-1)  becomes  indeterminate  since  both 
the  numerator  and  the  denominator  vanish.  The  behavior  of  the  solution  curve 
near  A  can,  nevertheless,  be  established  as  indicated  below.  Letting  z  =  1  +  J.' 
and  keeping  only  the  largest  terms  with  V  «1,  Equation  (n-l)becomes 


to  A  which  must  be  determined  by  numerical  integra- 
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d.i  [4  +  (v-l)(y-l)]  V  +  2 
dV  =  (n-l)V 


(n-io) 


with  the  singularity  now  at  /  =  0,  V  =  0.  Equation  (II -10)  is  linear  and 
can  be  solved  to  determine  the  behavior  of  the  solution  near  the  nodal  point  A. 
II  -  s.  when  V  =  V^then  in  the  cylindrical  case  with  v  =  2 

r2 

/  =  -  a  — 2  +  (y  +  3)V(  —  -  1)  (u-11) 

V  r 

r 


and 


lim 

V  ->0 


-  (7  +  3) 


i.e-,  the  solution  approaches  A  with  a  finite  slope.  In  the  spherical  (gse 
with  v  =  3 


V  V 
=  (y+l)V  In  —  -  a  — 

r  r 


(H-12) 


and 


lim 

V  -^>0 


d \X, 
aV 


_  rT1 


i.  e. ,  the  solution  approaches  A  with  infinite  slope. 

Numerically,  Equation  (II-l)  was  integrated  using  a  fourth  order 
Rur.ge-Kutta  Method.  The  integration  was  initiated  at  the  point  (Zg.Vg). 
which  is  known  once  the  C-J  conditions  have  been  determined.  Because 
of  the  singularity  at  A,  the  integration  was  only  extended  from  to  a 
minimum  value  of  0.005  for  V.  This  procedure  provided  an  adequate  repre¬ 
sentation  of  z(V)  even  near  the  singular  point  A.  Gome  typical  solution 
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curves  for  v  -  2,  3  are  shown  in  l-'igurc  (n  -1).  The  planar  case  m-1  is 
special  since  then  the  appropriate  solution  of  Equation  (n-1)  is  ihe  singular 
solution 

z  =  (1-V)2  (n-13) 

representing  a  plane  expansion  behind  the  detonation  front. 

As  indicated  above,  oncp  the  variations  of  z  with  V  is  determined,  z 

and  V  can  be  related  to  \  =  r/r  by  integration  of  Equation  (II -4).  Then  the 

s 

density  variation  can  be  found  from  Equation  (II-  5).  Since  the  flow  behind 
the  detonation  front  is  isentropic  the  pressure  ratio  (p/p2)  is  related  ti  the 
density  ratio  (p/p2)  by 


-E_  = 
p2 


while  the  temperature  ratio 


(T/Tg)  is  given  by 


(11-14) 


(r-)(7£) 


(11-15) 


Equations  (11-14)  and  (n-15)  are,  of  course,  based  on  the  assumption  that 
the  combustion  products  can  be  treated  as  a  perfect  gas  with  constant  specific 
heats. 


As  a  check  on  the  computational  technique  described  above>results  for 

/  (4) 

v  =  3,  p^  =  0,  =  y 2  =  5/3,  were  compared  with  the  results  of  Sedov  who 

also  considered  this  case,  and  exact  agreement  was  found.  As  a  typical 

example  pressure,  velocity  and  temperature  profiles  for  a  methane  air 
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detonation  are  shown  in  Figures  (II-2),  (E-3),  and  (II-4). 

The  impulse  functions  S( y y)  can  be  computed  once  the  variation 
of  p/p2  with  A  is  known.  From  the  definition  of  P  in  Equation  (n-3),  it 
follows  that 


P  = 


P 


Pl(t2) 


l_  p_  _ 


r2  P2 


P1  ,2 


— * L_  (P_) 
n  r2,2  Po 

pj  V  A  z 


(n-16) 


From  the  C-J  conditions  (Equation  (II-6)),  it  can  then  be  shown  that 


pic 


1  1 
= - r(i  +  — ) 

V1  >7 


fa -IT) 
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1  1 
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and  from  Equation  (25) 
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Figure  II-2.  Air-CH.  C  J  Detonation  Vari 
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As  indicated  in  the  main  body  of  this  report,  for  MAPP,  air 
and  methane-air  detonations  the  variations  of  y^,  y ^  and  are 

sufficiently  small  that  v )  essentially  depends  only  on  the  geo¬ 


metric  factor  v. 


APPENDIX  III 


CALCULATIONS  OF  CHEMICAL  EQUILIBRIUM  IN 
CHAPMAN- JOUGUET  CONDITIONS 

The  calculation  of  chemical  equilibrium  is  the  essential  feature  of 
the  Gordon  McBride  Program  (6).  This  calculation  is  carried  out  by 
iteratively  minimizing  the  Gibbs  or  Helmholtz  free  energy  as  described 
in  detail  in  Ref.  6.  The  program  computes  the  equilibrium  composition 
of  gas  mixtures  with  each  component  satisfying  the  perfect  gas  equation. 

With  the  equilibrium  calculation  as  a  base  the  program  can  carry 
out  the  following  calculations: 

(1)  Chemical  equilibrium  for  assigned  states  (T,  P) ,  (H,  P) ,  (S,  P) , 
(T,  V), (U,  V)  or  (S,V). 

(2)  Theoretical  rocket  performance  for  both  equilibrium  and 
frozen  compositions  during  expansion. 

(3)  Incident  and  reflected  shock  properties. 

(4)  Chapman -Jouguet  detonation  properties. 

Condensed  species  as  well  as  gaseous  species  can  be  considered. 

The  program  includes  thermodynamic  data  for  62  reactants  and 
421  reaction  species  in  the  form  of  coefficients  for  polynomials  fit  to 
the  data  by  the  method  of  least  squares.  The  data  is  taken  mainly  from 
the  JANAF  tables  (JANAF  Thermochemical  Tables.  Dow  Chemical  Co. , 
Midland,  Mich.,  Dec.  31,  1960  to  June  30,  1970.  Also  Ser.  A,  June 
1963;  Per.  B,  Jan.  1964;  Ser.  C,  April  1965;  Ser.  D,  Mar.  1966; 


The  program  input  specifies  the  oxidizer  and  fuel  composition 
enthalpy  and  density  in  the  case  of  the  C-J  option.  Certain  code  words 
specifying  which  option  is  to  be  used  must  also  be  specified. 

The  output  tabulates  the  properties  of  the  burned  and  unburned  gas 
and  the  detonation  parameters  p2/p1,  Tg/Tj,  p  / 'p  ,  and 

detonation  velocity  C.  The  program  also  provides  the  final  equilibrium 
composition  and  a  list  of  the  products  considered  in  the  calculation. 

A  typical  input  and  output  for  the  calculation  of  the  properties  of  a 
methane -air  detonation  is  reproduced  below.  Detailed  instructions  for 
use  of  the  program  and  sample  inputs  and  outputs  are  presented  in 
Ref.  6. 
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